G> . Standard Sources of Particle Production in Heavy Ion 

^ ' Collisions 

^-^ 

A. Capella 

O 

Q ; Laboratoire de Physique Theorique (CNRS - UMR 8627), 

I ■ Universite de Paris XI, Batiment 210, 91405 Orsay Cedex, France 

>. 

0^ ■ We describe particle production in the framework of an independent 

string model : the dual parton model. We show that an improved version 
of this model, containing a diquark breaking component, allows to describe 
the bulk of particle production and, in particular, baryon stopping and 
^y. \ most of the observed enhancement of strange baryons. Only for very rare 
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with final state interaction (comovers interaction). Recent data on event- 
by-event fluctuations in px are also described by the model. Predictions 
for RHIC and LHC arc presented and the effect of nuclear shadowing is 
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1. Introduction 



The enhanced production of strange particles (in particular of strange 
baryons) and the spectacular suppression of J/ijj in PbPb collisions are 
considered by many authors to be signals of Quark Gluon Plasma (QGP) 
production -or at least of the production of a system in thermal equilibrium. 
In this paper, we describe these phenomena in the framework of the dual 
parton model (DPM). We show that the model in its original form fails to 
describe the large amount of stopping measured in central SS and PbPb 
collisions -as do most independent string models. An improved version 
of the model containing a diquark breaking component is introduced. It 
allows to describe stopping without any extra free parameter. The same 
component is also responsible for most of the observed enhancement of 
strange baryons. However, the fi yield is underestimated by a factor five. 
Agreement with experiment is restored by introducing final state interaction 
(comovers interaction) with a very small cross-section of the order of 0.1 

(1) 



ARTICLERED PRINTED ON FEBRUARY 1, 2008 



mb. Comovers interaction is also needed in order to describe the anomalous 
suppression of J/ip in central PbPb collisions. 

Comovers interaction is a very non-trivial phenomenon mostly at a par- 
tonic level, which is not entirely understood. It turns out, however, that the 
cross-sections required to describe the data, both for strangeness enhance- 
ment and J/V' suppression are very small and the comovers interaction does 
not affect the bulk of particle production. In particular BB annihilation in 
the final state seems to be negligibly small. 

Recently, a lot of attention has been devoted to the study of event-by- 
event correlations -in particular in pT- The small value of this correlation 
in central PbPb collisions has been interpreted as a sign of thermalization. 
We show that this value is well reproduced in DPM. 

Predictions of the model for RHIC and LHC are also presented. Al- 
though minijets are important in order to determine pT distributions, in 
DPM they do not affect multiplicities. (For the latter, the short q-q strings 
in DPM play the same role as minijets). On the contrary, shadowing correc- 
tions are very important and reduce the values of dN/dy at mid-rapidities 
by a factor 2 at RHIC and by a factor 3 at LHC. This reduction is much 
stronger than in models where minijets are the dominant component. The 
reason being that in DPM shadowing corrections are present both for soft 
and semi-hard production. 

2. The model 

The dual parton model (DPM) is a dynamical model for low pT hadronic 
and nuclear interactions. It is based on the large-A^ expansion of non- 
perturbative QCD in the Veneziano limit -in which the ratio Nc/Nf is held 
fixed [1]. The dominant configuration consists in the production of two 
strings (of type qq-q in pp scattering, see Fig. 1). There are also more com- 
plicated configurations, corresponding to higher order terms in the large- A^ 
expansion, involving 4, 6, ... etc strings. These extra strings have sea quarks 
and antiquarks at their ends (Fig. 2). These configurations correspond to 
multiple inelastic scattering in an 5'-matrix approach. In pp collisions the 
contribution of each configuration is suppressed by a factor jV"^"'"'"^, where 
n is the number of strings -irrespective of the number of exchanged gluons 
and q-q loops, which do not change the topology of the graph. It is not 
known how to compute the numerical values of the various contributions 
from the QCD lagrangian. However, Veneziano has shown that there is a 
one-to-one correspondence between the various terms in the 1/A^ expansion 
and those in a multiple scattering theory (the number of strings being equal 
to twice the number of inelastic collisions). In view of that, we determine the 
contribution of each configuration to the total cross-section using a multiple 
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scattering model : generalized eikonal or perturbative reggeon calculus in 
hadron-hadron collisions and Glauber- Gribov model in collisions involving 
nuclei. 

For A — B collisions (for simplicity we consider here A = B) m. the ap- 
proximation of only two strings per nucleon-nucleon collisions, the rapidity 
distribution of secondaries is given by [2] 



J 11 J- 4 y^ 

ay 

(1) 
Here N{y) are the rapidity distributions of the individual strings, ua is the 
average number of wounded nucleons of A and n is the average number of 
collisions. Both Tia and n can be computed in the Glauber model. For 
instance for an average collision (i.e. integrated over impact parameter), 
one has 



n 



A^aNN/TAB « A^'^ . (2) 



Note that the total number of strings is 2n, i.e. two strings per inelastic 
nucleon-nucleon collision. 

The interpretation of (1) is obvious. With Tla struck nucleons, we have 
at our disposal Tla diquarks of projectile and target {qqAp and qqAj., re- 
spectively) and as many valence quarks. This accounts for the first term in 
(1). The remaining strings : 2(n — ua) have to be stretched by sea quarks 
and antiquarks, because the available valence constituents are all included 
in the first term ; this accounts for the second term of (1). Of course we 
should combine the valence and sea constituents of the projectile with those 
of the target in all possible ways. However, for linear quantities such as 
multiplicities, each ordering gives practically the same result. 

We can see from (1) and (2) that, if all strings would have the same 
plateau height (i.e. the same value of A^(0)), the plateau height in an average 
AA collision would increase like A^'^ . However, at present energies, the 
plateau height of the qs-Qs-, strings is smaller than that of qq-q ones, and the 
first term of (1) dominates. One obtains in this way the celebrated wounded 
nucleon model introduced by our Krakow hosts [3]. At higher energies the 
contribution of the sea strings becomes increasingly important. Therefore, 
in order to make predictions for RHIC and LHC we have to introduce the 
multistring configurations in each nucleon-nucleon collision. If their average 
number is 2K (this number can be computed in a generalized eikonal model : 
one gets ^ ~ 2 at y/s = 200 and K f» 3 at ^/s = 7 TeV) the total number 
of strings is 2Kn and Eq. (1) is changed into 
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^ (y) = flA \N'^'^^P''^^^T{y) + iV''Ap-99A:r(y) + {2K - 2)iV«--' 

dy 
+ {n-nA)2KN'i'-^' . (3) 

The hadronic spectra of the individual strings N{y) is obtained from 
a convolution of momentum distribution function and fragmentation func- 
tions. Both are assumed to be universal, i.e. the same in all hadronic and 
nuclear interactions. This makes the model very predictive, in particular 
regarding the energy and vl-dependences. Moreover, momentum distribu- 
tion functions and fragmentation functions are determined to a large extent 
from known Regge intercepts [2] [4]. Finally, hadrons produced in different 
strings are assumed to be uncorrelated (string independence). This is a 
simplicity assumption which does not follow from the large-A'^ expansion. 

So far, we have only considered configurations with an even number 
of strings. What about odd string configurations. A configuration with a 
single 3-3 string is possible in some cases, i.e. when it is possible to annihi- 
late a projectile quark with the corresponding antiquark in the target (Fig. 
3). However, in this case flavor quantum numbers are exchanged between 
projectile and target and this contribution decreases like l-^/s (secondary 
reggeon exchange). 

A configuration with three 3-3 strings is possible in pp annihilation (Fig. 
4). In this case, baryon number is exchanged between projectile and target, 
and the corresponding contribution is also expected to decrease in 1/ ^/s. In 
any case, it is known experimentally that this contribution is very small at 
high energy. 

3. Baryon stopping 

In DPM, the net baryon production takes place from diquark fragmen- 
tation. The produced baryon is fast in average. The data for pp collisions at 
SPS and for peripheral AA collisions can be reproduced in this way. How- 
ever, in the case of central SS and PbPb collisions, a huge baryon stopping 
has been observed by both the NA44 and NA49 collaborations. It cannot 
be reproduced in the model. Actually most string models in their original 
form fail to do so. 

The origin of this problem resides in the association of the net baryon 
production with the diquark. This is not necessarily the case. Indeed, let 
us consider again the three string graph for pp annihilation of Fig. 4. Here 
the valence quarks and antiquarks are found in their corresponding hemi- 
spheres and yet no baryon or antibaryon is present in the final state. This 
indicates that baryon number can be independent of valence quarks. It also 
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shows that it can be transferred over large rapidity distances and annihi- 
late with the corresponding antibaryon number -very much in the same 
way as quark and antiquark annihilate in the one string diagram of Fig. 
3. Assuming that (JannihUation ~ l/\/S) One obtains for the rapidity distri- 
bution of baryon number (when it is not associated with valence quarks) 
da/dy ~ exp(— l/2Ay). 

The above picture has its theoretical justification in the works of Dosch 
[5] and Rossi and Veneziano [6]. These authors have constructed a gauge 
invariant state vector of the baryon which leads to a picture of the baryon 
made out of three quarks bound together by three strings which join in a 
point called string junction (SJ). In this picture it is possible to transfer 
the SJ over large rapidity distances -leaving the valence quarks behind. 
In what follows, this component will be denoted diquark breaking (DB) 
component - while the conventional one will be denoted diquark preserving 
(DP) component. 

Recent work based on the above or related ideas can be found in Refs. 
[7]-[13]. However , in order to understand in this way the huge stopping 
observed in central AA collisions, one has to understand why the DB com- 
ponent described above is small in pp and quite large in central AA collisions. 
In Ref. [10], it was argued that the A dependence of the DB component 
is stronger than that of the DP one. In this way, a satisfactory description 
of the pp and central SS data was obtained -and predictions for central 
PbPb turned out to be correct. This was achieved by introducing one free 
parameter -which determines the ratio of the DB over DP contribution in 
pp. However, this approach is not entirely satisfactory since it requires a 
sort of fine tuning, namely the size of the DB component in pp has to be 
small enough not to spoil the agreement with the data, and large enough to 
describe the central AA data thanks to its larger A dependence. 

In recent works [12, 13], we have introduced a different implementation of 
the DB mechanism which avoids any fine tuning -as well as the necessity of 
an extra parameter. Our first remark is that, although the DB component 
is definitely present in the case of a single pp collision, it gives rise to a 
three-string configuration (Fig. 5). This is not the dominant, two-string, 
configuration in the laige-N expansion and, therefore, is expected to be 
small. Since its presence is not required to describe the pp data at SPS 
energies, we take it equal to zero for simplicity. (Actually any value small 
enough to be in agreement with the pp data can be introduced but it will 
have practically no effect on the results for central AB collisions). The 
second important remark [10] is that in the case of two inelastic collisions 
per nucleon (for instance an incoming proton scattering inelastically with 
two nucleons of a nuclear target), the above topological suppression does 
not occur. Indeed, in this case the dominant configuration has four strings 
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for both the DP and the DB components (see Figs. 2 and 6). Therefore a 
natural assumption is that, in this case, the two contributions have equal 
weights (1/2). The generalization to the case of n elastic collisions per 
nucleon is not so obvious. The assumption we have made [12, 13] (to be 
checked by comparing with experiment) is that there is an equal probability 
(1/n) for the SJ to follow any of the n collisions. In one of them, the SJ 
will join a valence diquark and hadronize in the conventional (DP) way. In 
all other cases, it will hadronize according to the DB mechanism. 

With these assumptions, the rapidity distribution of the net baryon 
{B — B) in A A collisions can be written as 



dy n 



dNop , ,_ _ JNdb 



riA^ h {n-UA) 



dy dy 



(4) 



Here n is the average number of binary collisions and Tia the average number 
of participants of A (for simplicity we consider only the case A = B). The 
average number of collisions per wounded (or participant) nucleon in njfiA- 
The probability of the DP component is thus fiA/n and that of the DB one 
1 — UA/n = (n — nA)/n- The extra factor ua in Eq. (4) ensures baryon 
number conservation (see below) . Note that in the case of a single collision 
per nucleon [n = n^), only the DP component is present. The exact form 
of dNos/dy is given in Refs. [12, 13]. Its main feature is the exp(— l/2Ay) 
factor, discussed above. In Eq. (4), dN£)p/dy is the conventional diquark 
fragmentation component. All relevant formulae are given in Ref. [14]. 
Both components are normalized to 2 (upon integration in y). In this way, 
the net baryon yield is 2nA as required by baryon number conservation. 

Note that in pp collisions at high energy there is also more than one 
collision per nucleon due to unitarity. Therefore, Eq. (4) gives a well defined 
prediction for the increase of stopping in hadronic collision at high energies. 
An obvious prediction is that the stopping will increase with the event 
multiplicity -large event multiplicity corresponding to a large number of 
interactions and/or strings. Preliminary data for the p-p asymmetry in DIS 
at HERA are in qualitative agreement with this prediction. 

The results for the net baryon yield B — B in central SS and PbPb 
collisions are shown in Fig. 7. The PbPb data are well reproduced. The SS 
ones (which show a larger stopping than in PbPb) are not so well described. 
Note, however, that the discrepancy between the data and the model pre- 
dictions without DB component has been substantially reduced. Note also 
that there is no free parameter here. 
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4. Strangeness enhancement 

We have argued in the previous section that in central A A colhsions, a 
large number of net baryons are produced at mid-rapidities and that they 
are dominantly made out of the SJ plus three sea quarks (see Fig. 6). It 
is then obvious than a large number of net A, H and fi (i.e. an increase of 
these yields per participant) will also take place. As a matter of fact, this is 
the only possibility to produce net O. The experimental value of the ratio 
Ct/Vi ~ 0.4 in central PbPb collisions at mid-rapidities is very much in favor 
of the above picture. Moreover, there will also be a substantial increase in 
the yield of K^ associated to the production of A's. 

In order to obtain the absolute yields of strange baryons one can use 
trivial quark counting arguments together with the ratio S = 2s/ {u + d) of 
strange over non-strange quarks in the sea. Details of the calculations, using 
values of S in the range 0.2 -^ 0.3, can be found in ref. [13]. Very similar 
results are obtained [12] by adopting a more phenomenological attitude, 
namely, fixing the A, H and fi yields at mid-rapidities from the pBe and 
pPb data [15]. Their values in central AA collisions are then obtained using 
Eq. (4). 

Before stating our results it is necessary to recall the mechanism of 
antibaryon production in DPM. It consists of two terms. One of them cor- 
responds to the usual diquark antidiquark pair production in the string 
breaking process. This term scales with the number of participants. The 
second one corresponds to the presence of diquark-antidiquark pairs in the 
nucleon sea. This term scales with n — ua- 

The corresponding formula (for A = B) is 

3 = ^A j—^ + in-nA) ^ . (5) 

ay dy dy 

For details see [12, 13] and references therein. 

Our results are shown in Fig. 8 (dashed lines). The p and A yields 
are well reproduced. The H's are slightly underestimated. However, the 
O's are too low by a factor 5. In an attempt to describe the fi yield we 
have introduced the final state interactions : vrA^ -^ KA, irN -^ KT,, 
ttA -^ KS, vrS -^ KE, and vrH -^ KQ, plus the corresponding reactions 
for the antiparticles. They are governed by the gain and loss differential 
equations [16] 

-^ = ^CFki Pk{x) pi{x) -^aik pi{x) pk{x) . (6) 

'^ ^ K/ k 

The first term in the r.h.s. of (6) describes the production of particle i 
resulting from the interaction of particles k and H. with space-time densities 
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p{x) and cross-sections ake (averaged over the momentum distributions of 
the interacting particles. The second term describes the loss of particle i 
due to its interaction with particle k. The initial densities are the ones 
obtained without final state interaction and the averaged cross-sections are 
taken to be the same for all processes. For details see [12, 13]. The data 
are reproduced with a value of the cross-section as small as 0.14 mb (full 
lines in Fig. 8). Note that we do not consider the inverse reactions required 
by detailed balance. These reactions give a negligibly small effect since 
PkPA ^ PttPn, etc. For the same reason we have neglected strange exchange 
reactions such as ttQ ^^ Kr.. Although the averaged cross-section can be 
larger than the value used above, this is overcompensated by the values of 
the involved densities. Since po, <C ps (by a factor 20 at initial interaction 
time), ttVl -^ Kr. is disfavored as compared to vrH -^ KQ. Likewise, since 
PR ^ Ptti ^■^ ~^ ^^ i^ disfavored as compared to ttH -^ KQ. 



Note that final state interaction is by no means a trivial effect. First, 
it represents a departure from the idea of independent strings. Second, 
and more important, a large contribution to the integrals (6) comes from 
interaction times of a few fermi, close to initial time where the system is in a 
dense pre-hadronic state. Actually, Brodsky and Muller [17] introduced the 
concept of comover interaction as a coalescence phenomenon at the partonic 
level, in order to describe the final state interaction. It is therefore clear 
that a lot of theoretical uncertainty is introduced in this way. The important 
result, however, is that the cross-sections required to describe the data are 
very small and do not affect the bulk of particle production. 



In Fig. 9 we show the ratio B/B for pPb and for four centrality bins 
in PbPb. All these ratios decrease significantly between pPb and central 
SS collisions and also between central SS and central PbPb. Although 
this decrease does not contradict thermal or QGP models (the increase in 
these ratios at fixed baryochemical potential can be overcompensated by an 
increase of the latter), it is not easy to obtain and most global fits in the 
framework of those models do not reproduce it. A more crucial test of the 
thermal and QGP models is the ratio of different types of antibaryons. An 
important result of our model is a ratio A/p significantly smaller than one. 
This is in agreement with recent preliminary data on the p yield from the 
NA49 collaboration [18] -together with published data on the A yield by 
the WA97 one [15]. This is in contradiction with the prediction of a A/p 
ratio significantly larger than one given by J. Rafelski [19]. 
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5. J /ijj suppression 

^ The decrease of the ratio J/ if) over DY when the centrahty increases 
was proposed by Matsui and Satz [20] as a test of a deconfining phase 
transition. Shortly afterwards this decrease was found in OU colhsions by 
the NA36 coUaboration. It has also been found in SU and PbPb -as well as 
in pA collisions (see C. Gerschel, these proceedings). The presence of J/tp 
suppression in pA collisions, indicates that another physical mechanism is 
at work. Most authors consider that it consists in the interaction of the 
pre-resonant cc pair with nucleons of the nucleus. A fit of all existing data 
gives a value of 6 -^ 7 mb for this absorptive cross-section [21]. The same 
mechanism allows to describe OU and SU collisions but fails to describe the 
PbPb data -which have a larger ("anomalous") suppression. The latter has 
been interpreted as a sign of a deconfining phase transition [21]. However, 
an alternative explanation has been proposed, based on the idea of comover 
interactions (tt + J/ip -^ D + D + X) of the same type introduced in the 
previous section to describe Vt enhancement. It has been shown in [22] 
that with a cross-section of 0.6 mb one obtains in this way a reasonable 
description of all the data. At a quantitative level, however, the comover 
picture tends to slightly overestimate the J/ip suppression both in central 
SU and in peripheral PbPb collisions and to underestimate it slightly in very 
central PbPb collisions. However, this disagreement is rather small (about 
2a). Moreover, there are large uncertainties both in the theory and in the 
experiment. In particular, recent data on J/ifj production in pA collisions 
by the E866 collaboration [23] lead to a smaller value of the absorptive 
cross-section. If these data were confirmed, the agreement of the comover 
picture would improve. Indeed, by reducing the value of the absorptive 
cross-section, there would be "more room" for comovers in SU collisions. 

The discussion above refers to the true J/^ over DY data. The NA50 
collaboration has also presented the so-called minimum bias (MB) analysis 
in PbPb collisions. This refers to the ratio 

DY \MBj,,p \DY J theory ' ^ ' 

Here MB is the single particle inclusive cross-section for charged particles 
for minimum bias events -i.e. without requiring the presence of a J/tp 
or a dimuon pair in the final state. The first factor in the r.h.s. of Eq. 
(7) has been measured experimentally, while the second one is calculated 
theoretically. The advantage of this procedure is that the statistic is very 



^ This subject has been discussed in great detail in the lectures by C. Gerschel. The 
discussion in this paragraph is very sketchy and assumes that the content of her 
lectures is known. 



10 ARTICLERED PRINTED ON FEBRUARY 1, 2008 



high (and the statistical errors very smah). However, systematic errors do 
not cancel here. 

The ratio obtained in this way shows no saturation at large Et- On the 
contrary, it continues to decrease steadily at the highest available values of 
Et- This feature cannot be reproduced in a comovers approach, at least 
in its present version. One has to note, however, that one needs some 
theoretical assumptions in order to compute the theoretical ratio MB/DY. 
In particular one assumes that the tail of the Et distributions of MB and 
DY are identical. Since both distributions show a very steep fall off at the 
tail, this assumption plays a very crucial role in the determination of that 
ratio at very large Et- 

The data [24] can be reproduced in a deconfining scenario [25]. One has 
to introduce two deconfining phase transitions, a first one for the x ^^^ '4'' 
and a second one for the direct J/^. However, in this approach, the < pT > 
of the J/-0 versus Et has a decrease at large Et which is not seen in the 
data. On the contrary, the comovers scenario gives a saturation at large Et 
both for the ratio J/ip over DY and for the < pT > of the J/ip versus Et 
(see last paper in [22]). 

In conclusion, the present data on J/ip suppression are very interesting 
but their interpretation is not yet established. Hopefully, the forthcoming 
RHIC data will allow to clarify the situation (see section 6). 

6. Shadowing corrections and predictions at RHIC and LHC 

Using Eq. (3) with K ~ 2 at RHIC and K ~ 3 at LHC we obtain for 
central collisions at 7 TeV [26] 

— (7/*~0)~2000 (y*~0)~7900. (8) 

dy dy 

These results are obtained without taking into account semi-hard collisions 
(minijets). As discussed in the introduction, the latter do not affect the 
multiplicities, since the average number of strings is constrained by unitarity. 
The fact that some of the q-q strings can be the result of a semi-hard gluon- 
gluon interaction, will affect the intrinsic pT of the string ends and, thus, 
the pt distribution of produced particles. However, average multiplicities 
are practically unchanged. 

As explained in [26], the values in (8) are upper limits. A reduction 
in these figures is expected from shadowing corrections. In hard processes, 
shadowing corrections in the nuclear structure functions are well known. In 
our approach, however, these corrections are present irrespective of whether 
the process is hard or soft. Moreover, they have the same physical origin 
and are governed by the same equations in both cases. 
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The physical origin of shadowing corrections can be traced to the dif- 
ference between Glauber model and Gribov field theory. The space-time 
picture of the interaction is very different in the two cases. Let us consider 
hadron- nucleus collisions. In the Glauber model, we have successive (billiard 
ball type of) collisions, while in Gribov theory we have "parallel" collisions 
of different projectile constituents with target nucleons. A key result of 
Gribov [27] is that the h — A amplitude can, nevertheless, be written as 
a sum of multiple-scattering diagrams with elastic intermediate states — 
which have the same expressions as in the Glauber model. However, in 
Gribov theory, there are extra multiple-scattering diagrams which contain 
all possible diffractive excitations of the incoming hadron as intermediate 
states. At present CERN energies, these extra diagrams lead to corrections 
to the Glauber formula of the order of 10 -^ 20 % in the total cross-sections. 
However, their contribution to dN/dy is much larger and leads to a reduction 
in the figures in (8) by about a factor 2 at RHIC and a factor 3 at LHC 
[28]. 

It is well known that the size of high mass excitations of the initial hadron 
is controlled by triple Pomeron couplings. It has been shown in [29] that 
the values of the triple reggeon couplings determined from soft diffraction, 
allow also to describe hard diffraction measured at HERA. It is also well 
known [30] that the latter determines the size of shadowing effects in the 
nuclear structure functions at low x. From the analysis in [29] it follows that 
at a scale Q^ ~ 1 GeV^ and x ~ 10~^ (the x-value relevant at RHIC) the 
shadowing in the Ph structure function leads to its reduction by a factor 
0.7. At X ~ 10~^ (relevant for LHC) the corresponding reduction is 0.6. 
Squaring these values (in order to take into account shadowing correction 
in both projectile and target in the case of PbPb collisions), we obtain the 
reductions by a factor 2 at RHIC and 3 at LHC as stated above. Similar 
results are obtained [28] considering that the process is soft and computing 
the modifications to dN/dy at y* ~ resulting from the extra terms in the 
Gribov theory -using the standard value of the triple Pomeron coupling. 

The above considerations show that the average virtualities relevant for 
the calculation of the shadowing effects in dN/dy at y* ~ at RHIC and 
LHC are of the order of 1 GeV^ . In models where the dominant contribution 
is semi-hard the relevant average virtualities, obtained from perturbative 
QCD, are higher and the shadowing effects significantly smaller. 

With the above values of dN/dy at y* ~ it is possible to compute 
the J/V' suppression at RHIC and LHC resulting from the two mechanisms 
discussed in section 5, namely nuclear shadowing and comover interaction. 
The suppression resulting from the former mechanism depends only on the 
absorptive cross-section which is expected to be the same at all energies. 
On the contrary, the density of comovers depends strongly on energy and, 
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as discussed above, on the shadowing effects. 

The ratio of J /ij) over DY in a very central PbPb coUision at SPS energies 
over the corresponding ratio in pp colhsions obtained in Ref. [22] is 0.23. 
(This takes into account both nuclear absorption and comovers interaction). 
Without nuclear shadowing, the corresponding ratios at RHIC and LHC 
are 0.03 and 10~^, respectively. When shadowing corrections are taken into 
account the corresponding ratios are significantly larger : 0.11 and 0.02. 

7. Event-by-event fluctuations in px 

Let us define the quantity [31] 

TV 

Z = J2z, (9) 

where A^ is the total number of particles in a single event and 

Zi = PT,- < PT > ■ (10) 

Here py. is the pT of particle i in the event and < • • • > denotes the average 
over all events. The correlation 6 is then defined as 



<N> 



V<z^ > (11) 



where < z^ > is determined by mixing particles from different events. In 
Ref. [31] a simplistic superposition model was considered in which 

< Z^ >AA ^ < Z"^ >NN ^2) 

< N >AA ~ <N >NN 

In this model (p is the same in NN and AA. Experimentally, it has been 
found [32] that (/) decreases by a factor 3-^4 from NN to central PbPb. 
This decrease has been interpreted in [31] as a sign of thermalization (see 
however Ref. [33]. 

It turns out that the decrease of (p observed experimentally is also repro- 
duced [34] in a Monte-Carlo formulation of the QGSM [4] (a model which is 
very close to DPM). The rapidity distribution, multiplicity distribution and 
Pt distribution measured experimentally [32] are also well reproduced by the 
model [34]. The results for cj) are [34] : 4>pp = 9.0 MeV, (ppbPb = 2.4 MeV at 
SPS energies and (j)pp = 76 MeV, (ppbPb = 79 MeV at RHIC. As we see the 
values of predicted at RHIC are much larger than at SPS and practically 
equal in pp and central PbPb. 
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8. Conclusions 

We have shown that the main properties of particle production in hadronic 
and nuclear collisions can be described in a dynamical string model, the 
DPM or QGSM. These models do not incorporate non-standard sources 
of particle production -such as the formation of a thermally equilibrated 
system or a quark gluon plasma. 

The huge stopping observed in central SS and PbPb collisions requires 
a modification of the model consisting in the introduction of a diquark brea- 
king component. The PbPb data are well reproduced in this way with no 
extra parameter and predictions for stopping in other systems are obtained. 
Strange baryon enhancement can also be described in the framework of this 
improved version of DPM. However, for rare processes such as $7 production 
(and to a much lesser extent for H production) one has to introduce final 
state interaction (comovers interaction). This mechanism is also required 
in order to describe J/ip suppression. Event-by-event fluctuations in pT are 
also well described by the model. 
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Figure Captions 
Fig. 1 : Two string diagram in pp. 

Fig. 2 : a) Four string diagram in pp. b) Four string diagram in pA. 

Fig. 3 : One string diagram in pp. 

Fig. 4 : Three string diagram in pp. 

Fig. 5 : Three string diagram for the diquark breaking component in pp. 
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Fig. 6 : Four string diagram for the diquark breaking component in pA. 

Fig. 7 : Rapidity distribution of the net baryon number B — B m. central 
SS and PbPb collisions. The full (dotted) line is the result with 
(without) the diquark breaking component. 

Fig. 8 : B + B yields at mid-rapidities for minimum bias pPb and central 
PbPb collisions at SPS energies. Full (dashed) lines are the results 
with (without) final state interaction. 

Fig. 9 : Same as Fig. 8 for the ratios B/B. 
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Figure 1 
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Figure 2 
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Figure 3 
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Figure 4 
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Figure 5 
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Figure 6 
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Figure 7 
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Figure 8 
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Figure 9 
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